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• Background and Aims There are many unresolved issues concerning the biochemistry of fructan biosynthesis. 
The aim of this paper is to address some of these by means of modelling mathematically the biochemical processes. 

• Methods A model has been constructed for the step-by-step synthesis of fructan polymers. This is run until a steady 
state is achieved for which a polymer distribution is predicted. It is shown how qualitatively different distributions can 
be obtained. 

• Key Results It is demonstrated how a set of experimental results on polymer distribution can by simulated 
by a simple parameter adjustments. 

• Conclusions Mathematical modelling of fructan biosynthesis can provide a useful tool for helping elucidate 
the details of the biosynthetic processes. 
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INTRODUCTION 

Lolium perenne ( perennial ryegrass) plants containing unusually 
high concentrations of the fructose polymer fructan (high sugar 
grasses, HSGs) were obtained in the 1970s by Pollock and Jones 
(1979). In recent years there has been growing interest in the pos- 
sible deployment of HSGs in grazed pastures (Turner era/. , 2006; 
Edwards et al, 2007; Parsons et ai, 2012) which may increase 
animal meat and milk production (Lee et ai, 2001; Miller 
et ai, 2001), decrease nitrogen (N) deposition to the soil and 
the generation by the soil of the greenhouse gas (GHG) nitrous 
oxide (Ellis et al., 2011), and have effects on methane 
(CH4) emanations from grazing animals (Ellis et ai, 2012). 
This could lead to substantial impacts on the environmental foot- 
print of pastoral agricultural systems. For these and other reasons, 
it is important to understand fructan biosynthesis and degrad- 
ation, including the dynamics and principal controls of the 
processes. 

Recently, we have developed improved LC-MS (liquid 
chromatography-mass spectrometry) methods to analyse fruc- 
tans in ryegrass (Harrison et al., 2009, 2011, 2012), and a 
series of observations prompted us to reinvestigate the 1-SST/ 
1-FFT hypothesis [1-SST is the sucrose: sucrose 1- 
fructosyltransferase enzyme, reaction (3) below; 1-FFT is the 
fructan: fructan 1-fructosyltransferase enzyme, reactions (4) 
and (5)]. Originally, the 1-SST/l-FFT model was proposed 
for the inulin-producing species Helianthus tuberosus 
(Edelman and Jefford, 1968), which has subsequently served 
as the conceptual basis for fructan biosynthesis in other 
higher plants. The 1-SST/l-FFT paradigm has been criticised, 



mainly because of discrepancies found between in vitro pro- 
duced fructan oligomer profiles and those present in planta 
(Cairns, 1993; Cairns et al, 2008, and references therein). 
Fructans in grasses are much more complex than the linear 
inulin fructans found in, for example, H. tuberosus due to the 
activity of additional fructosyltransferases (FTs) such as 
6G-FFT (Lasseur et al, 2006; Hisano et al, 2008) and 
6-SFT (Lasseur et al, 2011). These enzymes produce the 
fructan trimers 6G-kestose and 6-kestose which, when poly- 
merized by FTs, result in the synthesis of levans and the neo- 
series of inulins and levans. However, from a kinetic point of 
view all FTs catalyse the same reaction, i.e. the transfer of a 
fructose from a fructan oligomer to either sucrose or another 
fructan oligomer and can therefore be summarized into a 
single reaction equation. 

We have constructed a mathematical model of possible path- 
ways of fructan biosynthesis to test a range of assumptions and 
to compare predictions of the model with observations on 
fructan concentration and polymer length. Generally, the 
models are based on: (i) a constant supply of glucose (glc), (ii) 
the energetically neutral and reversible production of fructose 
(fru) from glucose [reaction (1)], (iii) the production of sucrose 
[reaction (2)], (iv) the production of 1-kestose (gf2) by 
the 1-SST reaction [reaction (3)], (v) the production of fructans 
of higher degrees of polymerization by successive transfer of 
fructose from kestose to the growing fructan polymer by FTs 
[reactions (4) and (5)], (vi) other possible transfers of fructose 
[reactions (6) and (7)] and (vii) the synthesis and transfer of 
small fructose polymers [reactions (8)]. 
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Mathematical models of reaction kinetics are used to 
make predictions of fructan production, including polymer 
length distributions, which are compared with our own measure- 
ments. This may lead to a better understanding of this poorly 
understood but important area of research. 

MATERIALS AND METHODS 

Plant materials and fructan extraction 

Loliumpereiine seeds from two lines ('Fennema' , 'PGl 13') were 
germinated and seedlings transferred to pots containing 
nutrient-rich potting mix. Plants were grown in controlled envir- 
onment chambers at two different temperature regimes with a 
14-h light and a 10-h dark period. Temperatures were set to 
either 10 °C constant or to 20 °C during the light and 10 °C 
during the dark period. Plants were regularly cut back (every 3 
weeks) and maintained as described previously (Rasmussen 
et a/., 2009). 

Blades (vegetative material above the ligule) were removed 
after a 3-week regrowth period, immediately frozen and 
ground in liquid nitrogen and freeze-dried. Ground plant 
powder (25 mg) was extracted twice in 1 mL 80 % ethanol and 
subsequently twice in 1 mL of water at 65 °C with constant 
shaking as described (Rasmussen et al, 2007). Extracts 
were centrifuged and supernatants combined for MS analysis 
of fructans. Combined supernatants were brought to dryness 
under vacuum, reconstituted into 1 mL of water, filtered 
through a 10-mm filter and transferred into high-performance 
liquid chromatography (HPLC) glass vials. 

Fructan analysis 

Filtrates (5 |jlL) were injected into the ultra-HPLC system and 
fructans separated on a Thermo Hypercarb column as described 
(Harrison et al. , 2009). MS analysis was performed using an LTQ 
ion trap mass spectrometer (Thermo Fisher Scientific, Waltham, 
MA, USA) with electrospray ionization in negative mode, and 
data collection over the mass range of 300-4000 allowing the 
collection of MS data for fructans up to DP (degree of polymer- 
ization) 49 (Harrison et al, 2009). Here, we report data corre- 
sponding to fructan oligomers in the DP range of DP 3-10 
only. As reported previously (Harrison et al, 2011) ions of 
fructan oligomers with a DP of 3-10 were predominantly 
singly charged, resulting in mass to charge ratios (m/z) for the 
deprotonated molecular ions [M-H] ~ of the individual DP oligo- 
mers of DP3 (=gf2) 503-3, DP4 (=gf3) 665-3, DP5 (=gf4) 
827-3, DP6 (=gf5) 989-3, DP7 (=gf6) 1151-3, DP8 (=gf7) 
1313-3, DP9(=gf8) 1475-3 and DP 10 (=gf9) 1637-3. 

FRUCTAN BIOSYNTHESIS: SCHEMES 
CONSIDERED AND SOME SIMULATIONS 

Various schemes are outlined, starting with the simplest. See the 
Appendix for details of the mathematics, parameterization of rate 
equations [e.g. eqns (A2), (A14)], and the rate : state equations 
defining the inputs and outputs for each biochemical species in 
the model (Thornley and France, 2007, pp. 21, 24). 

Throughout, it is assumed that glucose is provided at a constant 
concentration. Glucose (denoted by glc and also by g) is 



converted into fructose (denoted by fru and also by f ) reversibly 
according to eqns (A2) and (A14): 

glucose ^ fructose (1) 

Next glucose and fructose are converted irreversibly into sucrose 
(glucose-fructose, also denoted by gfl) by means of eqn (A3): 

glucose -t- fructose sucrose (gfl) (2) 

Energy (ATP) is generally required for this reaction (e.g. 
Thornley and Johnson, 2000, p. 301), although the possible 
modification of the kinetics by this requirement is not considered 
here. Any invertase or sucrase action, giving sucrose hydrolysis, 
is ignored. 

Scheme (a): the basic scheme 

Fructan biosynthesis in plants requires multiple, substrate- 
specific FTs. The currently most widely accepted hypothesis 
for the biosynthesis of fructans in plants is the 1-SST/l-FFT 
model proposed by Edelman and Jefford (1968). This hypothesis 
was originally designed to represent inulin biosynthesis in 
Jerusalem artichoke (Asteraceae), but is now generally applied 
to other plant systems as well. 

The first step in this reaction sequence is catalysed by the 
1-SST enzyme, which facilitates the transfer of a fructose unit 
from a donor sucrose molecule to the 1 position of the fructose 
on a sucrose acceptor molecule, forming a (32-1 glycosidic 
bond and resulting in a 1-kestose molecule (gf2) and glucose 
(Koops and Jonker, 1996; Van den Ende et al, 1996; Liischer 
etal, 2000a, b; Chalmers et al, 2003). This reaction is 

sucrose(gfl) -|- sucrose(gfl) — kestose(gf2) -|- glucose(glc) 

(3) 

Here, sucrose is denoted by gfl (glucose-fructose), kestose by 
gf2 (glucose-fructose-fructose) and the glucose monomer by 
glc. See eqn (A18) for the rate equation. 

The product of this reaction, kestose (gf2), can be used as 
donor and acceptor for a second fructosyl transfer reaction 
which is carried out by 1-FFT. This enzyme transfers a fructose 
unit from the donor 1-kestose molecule to the 1 -position of the 
fructose on an acceptor molecule, forming a (31-2 glycosidic 
bond, thus extending the fructose chain on the acceptor molecule 
by one fructose unit; glucose is always terminal. 1-FFT can only 
use fructans as donor molecules, but not sucrose, while sucrose, 
fructans and fructose itself [reactions (8)] can act as acceptor 
molecules (Jeong and Housley, 1992; Koops and Jonker, 1996; 
Liischer etal, 1993a, b; Van den Ende et al , 1 996). The reaction 
for kestose as the acceptor and donor of fructose is [eqn (A25)]: 

kestose (gf2) -t- kestose(gf2) — sucrose (gfl) -I- gf3 (4) 

where gf3 denotes the polymer glucose-fructose-fructose- 
fructose. This reaction can be generahzed [eqn(A26)] to 

kestose(gf2) -I- gfn sucrose(gfl) -I- gf(n -I- 1). (5) 

Here n is an integer {n = 2, 3, . . . , 9). gfn is a fructan with n 
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fructose units attached linearly to glucose. Here we consider 
polymers up to gflO. 

Forage grasses such as Lolium perenne (perennial ryegrass) 
accumulate a mixture of fructan types, namely the inulin and 
levan series, and the inulin and levan neoseries (e.g. Heldt, 
1997, p. 241). Here we just deal with the inulin series, as in (5). 
These series differ in the glycosidic bonds employed ((32-1 in 
inulins; (32-6 in levans), and also in the position of glucose in 
the fructan chain (terminal in the inulin and levan series; internal 
in the inulin and levan neoseries). The biosynthesis of these fruc- 
tans requires additional enzyme activities such as 6G-FFT 
(synthesizes fructans with internal glucose), and 6-SFT (synthe- 
sizes levans). The general assumption is that these FTs, like 
1-FFT, also transfer single fructose units and can only use fruc- 
tans (gfn; n > 2) as donors, but not sucrose. 

The basic scheme [Scheme (a)] comprises reactions (l)-(5). 
FTs use only gf2 (kestose) as a fructose (f ) donor; they transfer 
a single fructose molecule at a time, as in reaction (5) with 
n>2. Reaction (5) is applied for n = 2, 3, . . . , 9. The reactions 
proceed, of course, beyond n = 9, but we only programmed 
the problem as far as « = 9 as this suffices to describe the es- 
sentials of the problem and also our measured data. The time 
course is given in Fig. 1. Kestose (gf2) and higher polymers 
(gf3, gf4, . . . ) all overshoot to a decreasing extent before 
approaching the steady state. At the steady state (Fig. 2) each 
fructan polymer (gfn, n > 3) reaches the same concentration, 
which is half that of kestose (gf2). With glucose concentration 
glc, constant at 01 mol L~\ the steady-state concentrations 
are (mol L~'): fructose [fru] = 0-0333*; sucrose [gfl] = 0-09; 
kestose [gf2] = 0-023; and [gfn], n = 3, 4,..., 9 = 0-0112 
(Fig. 2). 

Scheme (b): extra fructose transfer donor added 

Modifying (a) above, it is now assumed that FTs can use both 
gf2 (kestose) [as in the basic scheme (a)] but also gf3 
(glucose-fructose-fructose-fmctose) as a fructose donor with 
transfer of a single fructose molecule. Thus, in addition to reac- 
tions (l)-(5), the reaction [eqn (A35)] 



0-10 



_ 0-08 
o 



1= 0-06 
o 

■£ 0-04 
o 

J 0-02 



glc (glucose) 



gfl (sucrose) 



/ 

/ 

I 

I 

I , .. 

/ / 



fru (fructose) 
gf'2 (kestose) 



2 3 
Time, f (d) 



4 



Fig . 1 . Scheme (a). Time course for scheme (a): the basic scheme, comprising 
reactions (l)-(5), with default parameters (Table 1). In scheme (a) transfer of a 
single fructan only occurs from kestose (gf2). Mathematical equations are 
given in Appendix A [eqns (A2), (A14), (A3), (A18), (A25) and (A26)]. gfn 
(« = 3, . . . , 9) denotes a fructan with ii fructose units attached to glucose; gf4- 
gf8 are not labelled due to lack of space. 
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Fig . 2 . Scheme (a): the steady state. Steady-state concentrations of fructans pro- 
duced by reactions (l)-(5) as in the basic scheme (a) with default parameteriza- 
tion (Table 1 ) are shown. Mathematical equations are given in Appendix A [eqns 
(A2), (A14), (A3), (A18), (A25) and (A26)]. In Scheme (a) transfer of a single 
fructose only occurs from kestose (gf2). glc denotes glucose; fru, fructose; gf«, 
a fructan with n fructose units attached to glucose; gfl is sucrose; gf2 is kestose. 



gf3 + gfn kestose(gf2) + gf{n -M), n = 3, 4, . . . , 9 (6) 



is included. Comparing Fig. 3 with Fig. 2, the steady-state con- 
centrations of individual gfn (n > 3) are now, relative to gf2 
which is higher, much lower. With glucose concentration 
([glc]) constant at 0-1 mol L~', the steady-state concentrations 
are (mol L"'): fructose [fru] = 0-0333*; sucrose [gfl] = 0-09; 
kestose [gf2] = 0-034; [gO] = 0-0068; and [gfn] (n = 4-9) = 
0-0062 (see Fig. 3). 

Scheme (c): effect of a fast reaction 

This scheme simulates the production of fructans as described 
under scheme (b) (Fig. 3), but where one of the reactions of reac- 
tion (6), namely that with n = 4: 



gf3 + gf4 kestose (gf2) + gf5 



(7) 



has its maximum velocity (v34_25) set to a high value, that is 
''34_25 = 10 [eqns (A35) with « = 4]. The result of this tenfold in- 
crease is a strong depression of the concentration of acceptor gf4 
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[reactions (l)-(5), Figs 1 and 2] is supplemented by the transfer of a single fruc- 
tose from gf3 (glucose-fructan-fructan-fructan) [reaction (6), eqn (A35)]. 
Steady-state concentrations are given, glc denotes glucose; fru, fructose; 
gfl, sucrose; gt2, kestose; gfn (;! = 3, 4, . . . , ), a fructan with n fructose units 
attached to glucose. Default parameters (Table 1) are modified with: 
V3«_2(„+i) = 1, n = 3, . . . , 9. 
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Fig. 4. Scheme(c):effectofafastreaction.Steady-stateconcentrationsoffractansareshown.ParameterizationisasinScheme(b)wherev3„_2(n+ = 1," = 3, ... ,9 
in reaction (6) [Fig. 3, reactions (l)-(6), eqns (A2), (A14), (A3), (A18), (A25), (A26) and (A35)]. However, now a liigh maximal velocity is assigned to reaction (7), 
namely V34 25 = 10 [eqns (A35) with « = 4]. glc denotes glucose; fru, fructose; gfi, sucrose; gf2, kestose; gf«, a fructan with n fructose units attached to glucose. 



(Fig. 4). That is, a gap is generated in the distribution of fructans 
according to degree of polymerization. With constant glucose 
([glc]) concentration of 0- 1 mol L~ \ the steady-state concentra- 
tions (mol L~^) of the other variables are: fructose [fru] = 
0-0333*; sucrose [gfi] = 0-90; kestose [gf2] = 0-037; [gf3] = 
0-0063; [gf4] = 0-0025; [gfn] = 0-0058, n = 5, . . . , 9 (Fig. 4). 



Scheme (d): fructose transfers to fructose polymers 

The occurrence of oligomeric carbohydrates containing ex- 
clusively fructose units (without any glucose) has been described 
in early studies on the carbohydrate composition of monocotyle- 
dons (see Archbold, 1940, and references therein). We therefore 
added an additional scheme (d), with putative transferase 
enzymes catalysing the following reactions: 

gf2(kestose)-|-f (fructose) gfl(sucrose) 

-|- f2(fructose — fructose) 
gf2(kestose) + f2(fructose — fructose) (8) 

gfi (sucrose) -I- f3 
gf2(kestose) -|- f3(fructose — fructose — fructose) gf5 

In the first two reactions, kestose (gf2) donates a single fructose, 
which can be accepted by fructose (f) or by f2 [eqns (A27)], 
yielding sucrose as a product. In the last reaction, an enzyme cat- 
alyses transfer of gf2 (kestose) to f3 [eqns (A28)], giving a gf5 
molecule with no other product. 

The basic reaction set of scheme (a) [reactions (l)-(5), eqns 
(A2), (A14), (A3), (A18), (A25) and (A26)] is applied as usual 
but one of the reactions, namely reaction (4) catalysed by 
1-FFT, is made zero by setting the velocity for this reaction to 
zero: V22 13 = 0 [eqn (A25)]. There is no production of gf3 or 
therefore gf4 [reaction (5) with n = 3]. 

At the steady state, concentrations of kestose (gf2) were very 
low and gf3 and gf4 were not produced at all, making a hole in 
the polymer distribution. Fructans with gfn (n = 5-9) and f2 
and f3 reached high concentrations, comparable to those of 



sucrose (gfi). With constant glucose concentration, [glc] of 
0-1 mol L~', steady-state concentrations of the other variables 
are (mol L"'): fructose [fru] = 0-031; sucrose [gfi] = 0-085; 
kestose [gf2] = 0-0087; [gf3] = 0; [gf4] = 0; [gf5-9] = 
0-031; [f2] = 0-031; [f3] = 0-031 (Fig. 5). 

Scheme (e): comparison with data 

As described above we grew ten genotypes from two 
L. peretine lines, the European cultivar 'Fennema' (F) and an ex- 
perimental breeding line 'PG113' (P) in two temperature regimes 
(one at constant 10 °C - 10/10, the other at 20 °C during the light 
and 10 °C during the dark period - 20/10). In blades (tissue 
above the ligule) harvested after 21 d of re-growth after defoli- 
ation 'PG113' accumulated significantly higher concentrations 
of total water-soluble carbohydrates (mg g~' d. wt) compared 
with 'Fennema' in both the 10/10 (means P = 322-2, 
F = 236-3; P < 0-0001) and the 20/10 (means P = 341-1, F = 
258-6; P< 0-0001) treatments (Rasmussen et al, 2009). We 
chose 21 d post-defoliation blade material as it had been shown 
previously that at this time the expression and activity of exohy- 
drolases is negligible and that fructan polymerization is preva- 
lent in these tissues (Morvan et al., 1997; Lasseur et al, 2007; 
Lee et al, 2010; Tamura er a/., 201 1). 

To determine the distribution of fructan oligomers with differ- 
ent DP, extracts of the above material were analysed by ion trap 
MS (Harrison etal, 2009). For most of the individual fructan oli- 
gomers, we detected several peaks corresponding to the miz of 
the individual fructan oligomers. Specifically, we detected 5 
peaks with the mIz of 503-3, two of which represent the non- 
fructan sucrosyl-galactosides raffinose and loliose (Amiard 
et al, 2003), the other three the fructan oligomers 1 -kestose, 
6-kestose and 6G-kestose (Harrison et al, 2012). The latter 
three are represented in Fig. 6 as the sum of fructan trimers 
[gf2 (DP3)]. Four peaks each with mIz ratios of 665-3, 827-3, 
989-3 and 1151-3 were detected, representing inulin and levan 
(neo)series fructan tetra-, penta-, hexa- and heptamers, respect- 
ively. The sums of peak intensities of each mIz are shown as 
gf3 (DP4), gf4 (DP5), gf5 (DP6) and gf6 (DP7) in Fig. 6. 
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Fig. 5 . Scheme (d): fructose transfers to fructose polymers allowed. Steady-state concentrations of fructans are shown. The basic scheme (a) is assumed [reactions 
(l)-(5)] but with reaction (4) (gf2 + gfl — > gfl + gf3) switched off, V22 13 = 0 [eqn (A25)]. Additional reactions (8) are added for transfer of fructose (f) from gf2 
(kestose) to f (fructose) and to f2 (fructose-fructose) [eqns (A27)] and also for the gf2 (kestose) -|- f3 — > gf5 polymerization [3rd of reactions (8), (A28)]. Default 
parameters (Table 1) are modified with i'2fi 10= 1. '^2n_\n ~ 1 and ^'2f3_5 = 1 [(8), (A27), eqns (A28)]. glc denotes glucose; fru, fructose; gfl, sucrose; gf2, 

kestose; gffi, a fructan with n fructose units attached to glucose; f2, difructose; f3, trifructose. 
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Three peaks with miz 1313-3 [sum = gf7 (DP8)], one peak with 
miz 1475-3 [gf8 (DP9)] and two peaks with mIz 1637-3 [sum = 
gf9 (DPIO)] were also detected (Fig. 6). 



As can be seen in Fig. 6, the pattern of fructan oligomer distri- 
bution in blades followed a comparable trend for both L. perenne 
lines and at both temperature regimes (Fig. 6). The highest rela- 
tive intensities of DP3- 10 fructans in blades were detected for 
tri- and tetrameric oligomers (DP3 and DP4), followed by 
approx. 4-fold lower intensities of miz signals corresponding to 
DPS -8 fructan oligomers. The lowest intensities were detected 
for DP9 oligomers with slightly higher concentrations of DPIO 
fructans. 

Parameters of our model have been adjusted to see if it is pos- 
sible to obtain agreement with our recent measurements. In fact, 
it is very easy to obtain such agreement, as can be seen when 
Fig. 7 (model predictions) is compared with observations in 
Fig. 6. In Fig. 7, minimal parameter changes have been made 
to the simplest model scheme with default parameters 
[(Table 1; scheme (a); Fig. 2; reactions (l)-(5), eqns (A2), 
(A14), (A3), (A18), (A25) and (A26)]. The velocity parameters 
alone were adjusted to give the results shown in Fig. 6, focusing 
on the relative amounts of fructan trimers (gf2, glucose-fructose- 
fructose) through to gf9, a glucose with nine fructoses attached to 
glucose, and ignoring the concentrations of glucose, fructose and 
sucrose. We do not believe our particular parameterization 
(Fig. 7) is unique and there are likely to be other sets of parameter 
values that would do equally well. 



DISCUSSION 

This paper presents a 'proof-of-concept' position, demonstrating 
that constructing and simulating specific biosynthetic schemes 
for fructan biosynthesis may significantly assist in understanding 
these processes. The mathematical analysis and computations 
are straightforward and, arguably, such approaches could be 
part of the standard armoury of techniques which is brought to 
bear on these problems. 

More particularly. Figs 2-5 illustrate, using simple assump- 
tions, the range of responses of polymer length distribution 
which can be obtained. Note also that the simulations presented 
here deal with fructan biosynthesis alone, without the possibly 
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conflicting and confounding effects from turning on fructan 
degradation processes (a topic that could be given a similar 
treatment). 

Figure 7 demonstrates that our scheme is able to 'explain' mea- 
sured data (presented in Fig. 6), although many more measure- 
ments would be needed to support or refute any detailed 



scheme proposed with its mathematical and numerical assump- 
tions. The only parameter change we introduced is a 5-fold in- 
crease in velocity of the reaction leading to the synthesis of gf8 
and a 2-fold increase of the velocity leading to gf9. These new 
parameters resulted in high concentrations of tri- and tetrameric 
fructan oligomers, relatively lower concentrations of penta-, 
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Fig. 7. Scheme (e): comparison with data (shown in Fig. 6). Scheme (a) is applied [reactions (l)-(5);eqns(A2), (A14), (A3), (A18), (A25)and(A26)]. Only non-zero 
default parameters are altered. Changes to the default parameters (Table 1) are: V22 13 = 2, V24 15 = 2, V25 le = 2, V26_i7 = 2, V27 ig = 2, V'28_i9 — 10. ^29 no = 4. 
Steady-state concentrations are: [gf2] = [gf3] = 0-016, [gf4] = [gf5] = [gf6] = [gf7] = 0 0079, [gf8] = 0 0016, [gf9] = 0-0039 mol L"'. gin denotes a fructan 

with n fructose units attached linearly to glucose. 



Table 1. State variables with initial values, variables, and parameters (concn. = concentration, mol = mole = 1 gram molecule, L = 
litre = 0-001 m^). Default values are given, corresponding to Figs 1 and 2 [Scheme (a), reactions (l)-(5)]. Relevant equation and 

reaction numbers are indicated 



State variable 



Definition 



Initial values (mol L ) 



(a) State variables 

[glc] 

[fru] 

[gfl] 

[gf2] 

[gfn], n = 3, 4, 5, 6, . . . , 9, 10 

[G] 

[f3] 



Glucose concn. (held constant) (AlO) 
Fructose concn. (A 17) 
Sucrose concn. (A24) 
Kestose concn. (A33) 

Fructan with n fructose units [(A38), (A44), (A51), (A58), . . . , (A63), (A67)] 
Fructose dimmer (A70) 
Fructose trimer (A73) 



0-1 

0 

0 

0 

0 

0 

0 



Parameter 



Definition 



Value and units 



(b) Parameters 

^fru_glci ^glc_fru 

"glcfru_gfl 

-'^gflgfl_glcgf2 

*^gf2gfn_^flgf(n+l) 



^gf2f3_gf5 
^gOgfn_gf2gf(n+l) 



Michaelis-Menten constants 
Fructose to/from glucose (1), (A2), (A14) 
Sucrose synthesis (2), (A3) 
Kestose synthesis (3), (A18) 
gf3 synthesis (4), (A25) 
(5), (A26), n = 3, 4, . . . , 9 
Synthesis of fructose dimer and trimer (8), (A27) 
Last of reactions (8), (A28) 
Reaction (6), (A35), « = 6, 7, . . . , 9 



0-1 molL"' 
0-01 (molL"')^ 



0-01 (molL"y 
0-01 (molL"')^ 
0-01 (molL"y 
0-01 (molL~y 
0-01 (molL"')^ 



glc ^glc_fru 
''glcfru_gfl 
^gflgfl_glcgf2 
Vgf2gG_gflga 
^gf2gfn_gflgf(n+l) 
^gf2fl_gflf2> l-'gf2C_gflfi 
Vgf2t3_gf5 
"gOgfnjQgfln+D 



Fructose to/from glucose (1), (A2), (A14) 

Sucrose synthesis (2), (A3) 

Kestose synthesis (3), (A18) 

gf3 synthesis (4), (A25) 

(5), (A26), « = 3, 4, . . . , 9 

Synthesis of fructose dimer and trimer (8), (A27) 

Last of reactions (8), (A28) 

Reactions (6), (A35), « = 3. 4, . . . , 9 
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hexa-, hepta- and octamers, a very low concentration of gf8 
(DP9) and slightly higher concentrations of gf9 (DP 10). One 
way of achieving this higher velocity in plants is the activity of 
an additional FT with a high affinity for higher DP fructans. 
Such FTs have been described for Echinops ritro (Van den 
Ende et al, 2006) and Phleum pratense (Tamura et al, 2009). 
To date no such high DP FT has been isolated from L. perenne 
and our model and data might indicate that the identification of 
such an enzyme would help to explain the observed differences 
of in vitro synthesized fructan profiles obtained by assaying 
known native L. perenne FTs compared with plant fructan pro- 
files (Cairns et a/., 2008). 

In view of the possible importance of the high-sugar pheno- 
type in grassland productivity, in mitigating GHG emissions 
and in possibly promoting carbon sequestration, we believe 
that a detailed mechanistic understanding given by analyses 
along the lines presented here may enable us to better understand 
the contributions of the high-sugar traits to grassland ecosystems. 
A worthwhile aim is to represent such characteristics realistically 
in grassland ecosystem models (e.g. Thornley, 1998). 

In conclusion, it has been demonstrated that mathematical 
analysis and simulation, based on well-established biochemical 
kinetics and calculus, is able to play a role in furthering our 
knowledge of fructan biochemistry, with its potential import- 
ance to the environmental and economic consequences of grass- 
land agriculture. 
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towards a product pool y with units of mol x per litre per day. 
Inputs to pool V from the x ^ y reaction are designated y 
(mol y d~ ). Where two substrates (x and y) are required 
giving two products (z and u), the output flux is designated 
Oxy zu ™d is calculated as an output from one of the substrate 
pools (say x) which requires a stoichiometrically related flux 
from the second substrate pool (y). The input flux I„._zu is in 
units of mol m L~ ^ d~ ^ or mol z L~ ^ d~ ' into the u and z pools. 
Parameter values are given where they are introduced and are 
listed in Table 1 for reference. 



Glucose ( glc) pool 

Synthesis of fructans is driven by maintaining a constant 
concentration of glucose [glc], namely 



[glc] = 01 mol L"'= Olg mol litre"'. 



(Al) 



Outputs. The first output flux from the glucose pool is to fructose 
[the forward reaction (1)] (mol glucose L~' d~') 

_ [glc] 

^'^-'™"'^'^-'""[glc] + /rgkj„ 

Vgic Jru — 1 niol glucose ' d^ ' , /^gic_fru =0-1 mol glucose L~ ' 

(A2) 



APPENDIX. MATHEMATICAL DETAILS AND 
PARAMETERIZATION OF THE MODEL 



The equations given below were programmed in ACSL 
(Advanced Continuous Simulation Language, Aegis Research, 
Huntsville, AL, USA; version 11-8-4), an ordinary differential 
equation solver, using fourth-order Runge-Kutta integration 
and a time step of 0-03125 (=1/32) d. The time unit of 
the model is days (d). There were no problems in model imple- 
mentation. 

We step through the biochemical state variables pool by pool, 
giving the outputs from, the inputs to, and the differential equa- 
tion of each pool. A constant volume of reactant, 1 litre (L), is 
assumed. A molar concentration is 1 gram molecule (mol) per 
litre (mol L~'), and a flux of 1 mol glucose L~' d~' is 1 gram 
molecule of glucose per litre per day. 

Fluxes are calculated as outputs O from a substrate pool, using 
a mechanistic equation, usually the Michaelis-Menten equation 
or something similar for bisubstrate reactions (Thornley and 
France, 2007, pp. 107-109, 113-114). If there is no substrate 
in a pool, there can be no output. Apart from the system input 
[eqn (A8)], inputs / are the consequence of outputs O from 
other pools, sometimes applying a simple stoichiometry, e.g. 
[eqns (A5), (A6)]. 

For a single substrate (x) single product (y) reaction, fluxes are 
calculated and designated as outputs ((9j-_^) from substrate pool x 



This is a traditional Michaelis-Menten (MM) equation. The 
parameters are: maximum velocity Vgi^ fru achieved at high 
values of [glc]; and an MM constant Kgi,._[^ giving the 
half-maximal- velocity glucose concentration [glc]. 

The second output flux is due to the two-substrate reaction 
(2) for sucrose synthesis (denoted by gfl, for glucose-fructose): 



[glc][fru] 



Cglcfru_gfl — Vglcfru_gfl . , „ ^ , 

Lglcj[ rruj -|- Agicfm.gfi 
Vgicfru_gfi = 1 mol glucose, fructose L~'d~\ 



(A3) 



glcfru. 



l\2 



an =0-01 (mol L-^ 



Again, the parameters are: a maximum velocity, v and a MM 
constant, K this time for a two-substrate reaction. This 
output flux is removed from the fructose pool [eqns (A 15)] as 
well as equally from the glucose pool and it is put into the 
sucrose (gfl) pool [eqn (A19)]. 

Total output from the glucose pool is (mol glucose L~' d~') 



Oo\c = Oglc_fi-u + O. 



glcfru _gfl 



Inputs. There is an input from the fructose 
tion(l): 

^fru_glc — ^fru_glc 



(A4) 



glucose reac- 



(A5) 
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nil 



This is calculated as an output from the fractose pool in eqn Fructose (fru) pool 
(A14). 

There is a second input to the glucose pool from the 1-SST 
reaction for kestose synthesis (3). This is calculated as (mol 
glucose d~^) 



Outputs. The first output is to glucose via reaction (1) in the 
reverse direction: 



4fl-glc - 2^gflgfl-gf2glc 



^fru_glc — ^fru_glc 



[fru] 



(A6) 



[ fru]+/i:fi-u_j,ic 



^fru_2lc 



1 mol fructose L 



C'gngf i_gf2gic is an output from the sucrose pool (gf 1 ) [eqn (Al 8)] 
below. 

Adding these two input fluxes gives a total input flux of 



^^fru.gic = 0-1 mol fructose L ^ . 

(A14) 



This is the mirror image of eqn (A2). The second output is to 
sucrose synthesis, Ogicftu gf i , in eqn (A3) above. The third 
/fiu_gic + 4fi_gic- (A7) output, Ogfan gfiG, is for the synthesis of the fructose-fructose 

dimer (f2) [first of reactions (8)], calculated as an output from 
The glucose concentration is (by assumption) held constant. This the kestose (gf2) pool below [first of eqns (A27)] . Total output is 



is achieved by adding an external (ext) input flux of 

4xt_glc = Ogle — (^fru.glc + 4fl-glc) (A8) 

Total input flux is 

Igic — ^fru_glc + ^fl-glc + ^ext_glc = Ogle (A9) 



Differential equation. The rate of change of glucose concentra- 
tion is 



d[glc] 
dr 



= /gic - Ogle, with [glc] 



Jo 



d[glc] 
dr 



d/, and 



(AlO) 



f = 0, [glc] = 0- 1 mol L' 



O, 



fru 



o, 



fru_glc 



o. 



glcfiu_gfl 



o 



gf2fl_gflf2 



(A15) 



Inputs. There is a single input to the pool from glucose, with 

^fru = 4lc-fru = OgieJiu (A16) 

where the output from the glucose pool is calculated in eqn (A2). 
Differential equation. This is 



d[fru] 



dt 



— hu — Of I 



u, with [fru] = 

Jo 



d[fru] 
dt 



dt, and 



(A17) 



? = 0, [fru] = OmolL" 



In view of eqn (A9), d[glc]/df is zero, and the glucose concen- Sucrose (gfl)pool 
tration remains at its initial value of 01 mol L~\ 

Notethatthetotalinputtothesystemis(molglucoseL-M-') Outputs. It is assumed that there is a single output from 

the sucrose pool giving kestose (gf2) synthesis as in reaction 
(All) (3) (invertase or sucrase activity, promoting sucrose hydrolysis, 
is ignored). This flux is 



ext_glc 



There are no outputs from the system so the rate of change of the 
whole system variable {Sys mol glucose L~' d~') is 



dt ~ ' 



ext_glc 



(A12) 



A useful check on the formulation and programming of the 
problem is also to calculate the rate of change of the whole 
system by summing its parts: 



dS, d[glc] d[fru] d[gfl] 

(check) = — h 2— 

dt dt dt dt 

, ^d[gf8] , ^Qd[gf9] , 



dt 



dt 



d[gf2] 

dr 
d[f2] 
dt 



, d[f3] 
dt 



Equations (A12) and (A13) should be identically equal at 
all times. 



Ogfl — Ogfigfi_glegf2 — Vgfigfi_glegf2 



[gfl][gfl] 



[gfl][gfl]+^gflgfl_glcg 



Vgflgfl_glcgf2 



= 1 mol sucrose L d 



K2 



U-1 



%igfi_gicgf2 = 0-01 (mol L ) 



(A18) 



Inputs. The principal input to the pool is due to sucrose synthesis 
from glucose and fructose, reaction (2), with [eqn (A3)] 



^glcfru_gfl — Oglefru_gfl 



(A19) 



(A13) A second input is from the kestose pool, reaction (4), given by 



4f2gf2-gflgt3 = 2^gf2gf2_gflgf3 



(A20) 
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This output from the kestose pool is calculated in eqn (A25). synthesis of fractose dimers and trimers [reactions (8)] 
There are also many similar inputs, arising from the outputs of 
reactions (5) with n = 3, 4, . . . , 9; 



^gf2gfn_gflgf(n+l) 



o. 



gf2gfn_gflgf(n+l) 



(A21) 



These outputs from the kestose pool are calculated in 
eqn (A26). 

There are two more inputs, arising from the synthesis of 
the fructose dimer f2, and the fructose trimer f3, shown in reac- 
tions (8): 

4f2fl-gflf2 = Ogf2fl_gflf2 and /gf2f2_gflf3 = Ogf2G_gflD (A22) 

The outputs on the right of these equations are again calculated 
under the kestose pool, eqns (A27). 
Total input to the sucrose pool is 



O 



'gf2fl_gflf2 — I^gf2fl_gflf2 



[gf2][fru] 
[gf2][fru] + /i:gf2fi_gfif2 

Vgf2fi_gfit2 = 0 mol(gf2, fru)L"',/i:gf2fi_gfif2 = 0-01 (molL"')^ 

[gf2][f2] 



<3gf2f2_gflf3 = Vgf2f2_gflf3 



[gf2][f2] + /:gf2i2_gflf3 



Vgf2f2_gfio = 0 mol(gf2, f2)L-' , K^f2_inn = 0-01 (molL-i)^ 

(A27) 



Note that these reactions are turned off in default giving just the 
basic reaction scheme [scheme (a), reactions (l)-(5)]. 

There is one further potential output to be considered, that of 
the third of reactions (8). This is 



^gfi 



■ ■'glcfru.gfl + 4t2gf2_gflgf3 + 4f2gf3-gflgf4 

+ \- 4f2gf8_gflgf9 + 4f2gf9-gflgfl0 

+ 4f2fl_gflf2 + 4f2f2-gflf3 



(A23) 



:f2t3_gf5 = Vgf2f3_gf5 



[gf2][f3] 



[gf2][f3] + /:gp.f3_gf5 



Vgf2f3.gf5 = 0 mol L-'(gf2, f3) d"', /:gf2f3_gf5 = 0-01 (mol ir' f 

(A28) 



Differential equation. This is 
d[gfl] 



df 



/gfi -0^1, with [gfl] 
f = 0, [gfl] = OmolL^^ 



. r d[gfl]^ 
Jo df 



dr, and 
df (A24) 



Kestose pool ( gf2 ) 

Outputs. There is an output [reaction (4)] giving gf3 synthesis: 



<3gf2gf2gf lgf3 = I'gf2gf2.gflgf3 



[gf2][gf2] 



[gf2][gf2]+/:gf2gf2.gflgD 



Total outputs from the gf2 (kestose) pool are 

Ogf2 = Ogf2gf2_gflgf3 + Ogf2gf3_gfigf4 + •■■ - 
+ <?gf2gf9_gflgfl0 + C'gf2fl_gflf2 
+ <3gf2f2_gflf3 + Ogf2f3_gf5 



o, 



gf2gf8_gflgf9 



(A29) 



Inputs. The principal input to the kestose pool is via the 1-SST 
reaction from sucrose (3) with [see eqn (A 18)] 



4flgfl-glcgf2 = 2'^g"gfl-S''=g*^ 



(A30) 



Vgf2gf2.gfigf3 = lmolkestoseL d ,A'gf2gf2.gfigf3 = 0-01 (mol L ) 

( A25) There is also a series of inputs from reaction (6) where gf3 can act 
as a fructose donor [this is outside the basic scheme of reactions 
(l)-(5)]: 



There are a series of similar outputs from reaction (5) 
(n = 3, 4, . . . , 9): 



Cgf2gfn_gflgf(n+1) — Vgf2gfn_gflgf(n+1) 



[gf2][gfn] 



4f3gf3-gf2gf4 - 2<^gf3gf3-gf2gf4, /gf3gfn_gf2gf (n+1) 



[gf2][gfn]+/s:gf2gfn_gflgf 



Vgf2gfn_gfigf(n+i) = 1 mol (gf2, gfn) L M 

^gf2gfn_gflgf(n+l) = 0-01 (mol L~')^ 



(n+1) 



(A26) 



,f3gfn_gf2gf(n+l) 



.n = 4, 5,6, ...,9 (A31) 



The right side outputs are calculated under the gf3 pool below 
[eqn (A35)]. Total input is 



There can be outputs from kestose of a fructose monomer to a 
fructose monomer or to a fructose dimer, accomplishing the 



'gf2 = -'gflgfl_glcgf2 + ^gf3gf3_gf2gf4 + 'gf3gf5_gf2gf6 ' 
+ /gDgf8_gf2gf9 + 4f3gf9_gf2gflO 



(A32) 
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Differential equation. This is Total output from the gf4 pool is 



= Lf2 - 0,f2, With [gf2] = f dt, and '^gf4 = Ogf2gf4-gflgf5 + Cgf3gf4_gf2gf5 (A41) 

' Jo dr (A33) 

f = 0, [gf2] = OmolL^^ 



Inputs. There are two inputs [with eqns (A26) with n = 3 and eqn 
(A35) withn = 3]: 

Sf^P°°^ /gf2gf3_gflgf4 = Ogf2gf3_gflgf4, 4f3gf3-gf2gf4 = 2'^Sf3gf3-gf2gf4 

Outputs. There is an output of reaction type (5) with n = 3 giving /■ aai\ 

gt4 synthesis: 

Ogf2gO_gfigf4 (A34) Total input is 

This is given by eqn (A26) with « = 3. There is a series of outputs ■'gf^ — 4f2gf3-gfigf4 + 4f3gf3_gf2gf4 (A43) 

from reaction (6) (n = 3, 4, . . . , 9): 

Ogf3gf„_gf2gf (n+i, = VgOgf„_gf2gf („+i) . ^^ir,Jif/^^"^ Differential equation. This is 

Igr^J lginj+ Agf3gfn_gf2gf (n+ 1) 

VgOg.W(n.i) = 0 mol (gf3, gfn) L"'d-, d[g^ ^ ^^^^ _ ^^^^^ ^.^^ ^^^^^ ^ f' d[g^^^^ 

W-gf2gf(n+i) = 0.01(molL-')2 Jo dt (A44) 

(A35) f = 0, [gf4] = 0molL-' 

Total output from the gf3 pool is 

<?gf3 — C'gGgf3_gflaf4 + C'gf3gf3_gGgf4 + C'af3gf4 gf2gf5--- . , „ gf^ POol 

(Ajd) 

+ Ogf3gf8_gf2gf9 + Ogf3gf9_gf2gfio Outputs. There is an output of reaction type (5) and eqn 

(A26) with n = 5 to gf6 synthesis; 

Inputs. The only input is via the 1-FFT reaction from kestose (4) ^gf2gf5-gfigf6 (A45) 

with [see eqn (A25)] 

There is a second output from reaction (6) and eqn (A35) with 
1 n = 5: 

4*3 = 4f2gf2_gflgf3 = 2 '^gf2gf2-gflgf3 (A37) 

Ogf3gf5_gf2gf6 (A46) 



Differential equation. This is Total output from the gf5 pool is 

d[gf3] r ^ , [' d[gf3] 

(A38) 



- /gf3 - (3gf3 , with [gf3] = J ^ dr, and ^ ^ ^ Ogfs = Cigf2gf5_gfigf6 + C)gOgf5_gf2gf6 (A47) 



[gf3] Omol L Inputs. The first two inputs are [reaction (5), eqns (A26) with 

n = 4; and reaction (6), eqns (A35) with n = 4]: 

4f2gf4_gflgf5 = Ogf2gf4_gflgf5, 4f3gf4-gf2gf5 = <?gf3gf4_gf2gf5 

Sf4pool (^4g) 
Outputs. There is an output of reaction type (5) and eqn (A26) 

with n = 4 giving gf5 synthesis: There is a third input from the last of reactions (8) [eqn (A28)] 

Cgf2gf4_gflgf5 (A39) 4f2f3_gf5 = Cgf2f3_gf5 (A49) 

There is a second output from reaction (6) and eqn (A35) Total input is 
with n = 4: 

Ogf3gf4_gf2gf5 (A40) 4f5 = 4f2gf4_gflgf5 + 4f3gf4-gf2gf5 + /gf2t3_gf5 (A50) 
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Differential equation. This is 

= /gf5 - Ogf5, with [gf5] = — — df, and 

Jo 



dr 



dr 



(A51) 



f = 0, [gf5] = 0molL- 



gf6pool 

Outputs. There is an output of reaction type (5) [eqn (A26) with 
n = 6] to gf7 synthesis: 



O 



gf2gf6_gflgf7 



(A52) 



There is a second output from reaction (6) and eqn (A35) with 
n = 6: 



Ogf3gf6_gGgf7 

Total output from the gf6 pool is 

Ogf6 = Ogt2gf6_gflgf7 + <3gf3gf6_gf2gf7 



(A54) 



Inputs. There are two inputs. The first is [reaction (5); eqns (A26) 
withn = 5]: 



4f2gf5-gflgf6 = C'gf2gf5_gflgf6 

The second is [reaction (6) and eqn (A35) with n = 5]: 

4f3gf5-gf2gf6 = C>gf3gf5_gf2gf6 

Total input to the gf6 pool is 

Igf6 — /gf2gf5_gflgf6 + /gf3gf5_gf2gf6 



(A55) 



(A56) 



(A57) 



Differential equation. This is 
d[gf6] 



dt 



4f6 - Ogf6, with [gf6] 



Jo 



d[gf6] 
dt 



dt, and 



(A58) 



f = 0, [gf6] = 0molL" 



gf7 and gfS pools 

The equations are the similar to those for the gf6 pool [eqns 
(A52) to (A58) ]: for the gf7 pool, change 7 to 8, then 6 to 7 
and 5 to 6; for the gf8 pool, change 7 to 9, 6 to 8 and 5 to 7. 



gf9pool 

Outputs. There are two outputs: 

<5gf2gf9_gflgfl0> Cgf3gf9_gf2gflO 



(A59) 



The first is an output of reaction type (5) with n = 9 synthesizing 
gflO [eqn (A26) ]; the second is from reaction (6) and eqn (A35) 
with n = 9. Total output from the gf9 pool is 



Ogfi) — OgOgff.gflgflO + <3gf3gf9_gf2gfl0 



(A60) 



Inputs. The two inputs are [reaction (5), eqns (A26) with « = 8; 
and reaction (6), eqns (A35) with « = 8]: 



4f2gf8-gflgf9 = <3gf2gfS_gflgf9> 4f3gf8-gf2gf9 
= Cgf3gf8_gf2gf9 



(A53) Total input is 



^gf9 = 'gf2gf8_gflgf9 + 'gf3gf8_gf2gf9 



(A61) 



(A62) 



Differential equation. This is 
d[gf9] 



df 



/gf9-Ogff,with[gf9] 
f = 0, [gf9] = 0molL-' 



d[gf9] 
, dr 



dt, and 



(A63) 



The gf 9 pool is (in the present formulation) the last pool to reach a 
steady state (Fig. 1). Therefore, the proportional or relative 
growth rate of gf9, Rgfg (d~ is calculated to assess the accuracy 
of the final state. This is 



if([gf9] > 0) then Rgfg = 



1 d[gf9] 
[gf9] df 



endif 



(A64) 



gflO pool 

This is the last pool included in our analysis. It has inputs of 
[eqns (A26) and (A3 5) both with n = 9] 

4f2gf9_gflgfl0 = <3gf2gf9_gflgfl0, 4f3gf9_gGgflO 

= <?gf3gf9_gf2gflO (A65) 

giving a total input of 

tgfU) = 4f2gf9_gflgfl0 + 4f3gf9_gf2gfl0 (A66) 

There are no outputs so that gflO increases linearly in the steady 
state, when all the pools except gflO are constant [see eqn (A8) 
for the total system input, /eAt_gic]- Thus, in the steady state, the 
equation 



^ext_glc 



dt 



;flO 



(A67) 
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must be satisfied. This equation provides another useful check on f3 pool (fructose trimer) 



the accuracy of problem formulation and programming. 



f2 pool (fructose dimer) 

Outputs. The sole output is with the second of reactions (8) given 
in eqns (A27). Total output from the f2 pool is 



On — Cgf2f2_gflf3 



(A68) 



Inputs. There is one input, from the first of reactions (8) and the 
first of eqns (A27): 



'gf2fl_gflf2 



o 



gf2fl_gflf2 



(A69) 



Outputs. The sole output is with the last of reactions (8) given in 
eqns (A28). Total output from the f3 pool is 



On = O, 



gf2f3_gf5 



(A71) 



Inputs. There is one input, from the second of reactions (8) and the 
third of eqns (A27): 



gf2f2_gflf3 



= o 



gf2f2_gflf3 



(A72) 



Differential equation. This is 
d[f2] 



At 



In - On. with [f2] = | ^^At, and 



f 

JO 



At 



(A70) 



f = 0, [f2] = 0molL- 



Differential equation. This is 
d[f3] 



At 



= /f3 - (3f3, with [f3] = 



f 

Jo 



d[f3] 
At 



At, and 



(A73) 



f = 0, [f3] = 0molL- 



